Commensurate-Incommensurate Magnetic Phase Transition in Magnetoelectric Single 

Crystal LiNiP04 
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Neutron scattering studies of single-crystal LiNiP04 reveal a spontaneous first-order 
commensurate-incommensurate magnetic phase transition. Short- and long-range incommensurate 
phases are intermediate between the high temperature paramagnetic and the low temperature an- 
tiferromagnetic phases. The modulated structure has a predominant antiferromagnetic component, 
giving rise to satellite peaks in the vicinity of the fundamental antiferromagnetic Bragg reflection, 
and a ferromagnetic component giving rise to peaks at small momentum-transfers around the ori- 
gin at (0, ±Q,0). The wavelength of the modulated magnetic structure varies continuously with 
temperature. It is argued that the incommensurate short- and long-range phases are due to spin- 
dimensionality crossover from a continuous to the discrete Ising state. These observations explain 
the anomalous first-order transition seen in the magnetoelectric effect of this system. 

PACS numbers: 75.25.-(-z, 75.50.Ee, 78.20.Ls 



Common magnetic systems with simple colinear long- 
range-order ground state can melt into the paramagnetic 
(disordered) state directly, usually via a second order 
phase transition, or through a series of intermediate spa- 
tiallv modulated phases before losing all correlations P, 
0, 13 ■ The indirect melting through modulated phases 
indicates the presence of competing interactions of next 
nearest- neighbors, anisotropics in the spin Hamiltonian, 
and/or topological frustrations JJ. There has been a con- 
tinuous interest in the spontaneous and magnetic-field in- 
duced commensurate-incommensurate magnetic (C-IC) 
transition over the years^ |^ 0. For instance, in the 
semimetallic europium tri- arsenide (EuAss), the ground 
state of the system is commensurate and as temperature 
increases the system undergoes a C-IC transition In 
copper metaborate, on the other hand, the ground state 
is incommensurate and undergoes a continuous phase 
transition to a non-collinear commensurate antiferromag- 
netic state 0- It has also been demonstrated that the 
C-IC transition can be induced by the application of an 
external magnetic field0,|3. 

Here, we report a novel C-IC magnetic phase transi- 
tion in the weakly coupled antiferromagnetic planes of 
LiNiP04 (S=l, Ni^+), its characteristics resemble IC 
structural phase transition s [a. LiNiP04 is an antifer- 
romagnetic (AF) insulator [3, [ifl which belongs to the 
olivine family of lithium orthophosphates LiMP04 (M = 
Mn, Fe, Co, and Ni)[l3]; space group is Pnma 12]. Neu- 
tron scattering studies demonstrated that LiMP04 (M 
=Ni, Co) exhibit properties between the two-dimensional 
(2D) and three-dimensional (3D) with an interlayer cou- 
pling that is relatively stronger than the coupling found 
in the cuprates, for instancep^ IT^ . These insulators 
also exhibit a strong linear magnetoelectric (ME) effect, 
with the observed ME tensor components, Oixy,Oiyx, for 
LiCoP04 and, oixz^ctzxi for LiNiP04, in agreement with 
the antiferromagnetic point groups mmm' and mm'm. 
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FIG. 1: Magnetoelectric coefficients of LiNiP04 versus tem- 
perature measured by the dynamic technique. The ME coef- 
ficients, axz, and a^x, were measured under 5 kOe magnetic 
field along the c-axis and the a-axis, respectively ^^l- 



respectively, but with some anomalies^^ IT^ . In 
particular, the ME effect measurements of LiNiP04 as 
a function of temperature reveal a first-order AF tran- 
sition, and an unusual decrease of the ME coefficient 
at temperatures below a maximum close to Tn^- 
By contrast, the isostructural LiCoP04, LiFeP04, and 
LiMnP04 exhibit continuous change of the ME coeffi- 
cients, indicative of second-order transitions 15]. Mag- 
netic susceptibility studies of poly crystalline LiNiP04 
showed a significant deviation from the Curie- Weiss law 
in a temperature range much higher than T^r, and neu- 
tron scattering from the same polycrystalline sample gave 
rise to diffuse scattering at the nominal position of the 
AF Bragg reflection up to T « 2T/v 13] . Recent magnetic 
susceptibility measurements of single crystal LiNiP04 
showed two features, one at T^r = 20.8 K and one at 
Ti = 21.8 K associated with an AF transition and an 
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intermediate-phase, speculated to be IC, respectively 0. 

An irregular shaped single crystal (0.396 grams in size; 
lattice constants at RT: a = 10.0317 A, 6 = 5.8539 
A,c = 4.6768 A), synthesized by a flux method described 
elsewhere |23|, was used for the neutron scattering stud- 
ies. Neutron scattering measurements were carried out 
on the HBlA triple-axis spectrometer at the High Flux 
Isotope Reactor (HFIR) at Oak Ridge National Labo- 
ratory. A monochromatic neutron beam of wavelength 
A = 2.368 A (14.7 meV, ko = 2tt/X = 2.653A-1) 
was selected by a double monochromator system, us- 
ing the (0, 0, 2) Bragg reflection of highly oriented py- 
rolytic graphite (HOPG) crystals. The A/2 component 
in the beam was removed (to better than 1.3 parts in 
10^) by a set of HOPG crystals situated between the two 
monochromating crystals. The coUimating configuration 
40', 40', iSampZe, 34', 68' was used throughout the exper- 
iments. HOPG was also used as the analyzer crystal. 
Temperature measurements and control were achieved by 
a Conductus LTC-20 temperature controller using Lake 
Shore silicon-diode temperature sensors (standard curve 
10). ME effect measurements on thin polished plates 
with evaporated gold layers as electrodes were performed 
with the dynamic technique and quasistatic technique for 
calibration, as described elsewhere [T^. Figure shows 
a strong ME effect with an abrupt transition, with an 
anomalous temperature dependence of the ME coeffi- 
cients, azx and a^z, with a maximum close to Tn, in 
agreement with previous results fislll^ . 

Relatively wide longitudinal scans along the (0,K,0) 
direction at selected temperatures show the transition 
from the AF to the paramagnetic state proceeds through 
an infinite series of modulated structures. At low tem- 
peratures (below T « 19 K), a single Bragg reflection 
at T^p = (010), due to the colinear AF orderingT^, is 
observed, as shown in Fig.|5]for T = 9.81 K. At temper- 
atures higher than « 19 K and lower than T ~ 20.80 K, 
this Bragg reflection is superimposed on a very weak dif- 
fuse scattering in the form of a broad Lorentzian-shaped 
peak, also centered at (0,1,0). This broad peak (width 
0.235 A-i), due to short -range inplane coherence-lengths 
(^ « 27 A), is likely related to the lamellated domains 
recently observed by magnetic second-harmonic genera- 
tion topography on a thin plate of LiNiP04[^, and with 
the maximum observed in the ME effect (see Fig. At 
T = 20.80 K, two extra satellite reflections at (0,1±Q,0) 
appear, signaling a transition from the simple colinear 
AF phase to the IC magnetic phase. The transition oc- 
curs within 0.005 K of Tn (the resolution in temperature). 
As shown, at this temperature commensurate and incom- 
mensurate phases coexist, as is typical of first-order phase 
transitions and consistent with the sudden disappearance 
of the linear ME effect (former observations show that the 
linear ME-efFect cancels out in magnetically incommen- 
surate structures; e.g. BiFe03'22i, BaMnF4 jHI^). As 
the temperature is raised, the IC structure predominates 
with correlation lengths comparable to those of the long 
range AF ground state. In fact, long-range incommen- 
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FIG. 2: A) Longitudinal scans along the (0,K,0) direction 
showing a single sharp (resolution limited) peak at the (010) 
position (r.I.u is reciprocal lattice unit, in this case, b* = 
1.0783 units; the intensities of scans are shifted by two 
decades each for clarity). A single Bragg reflection, due to 
the AF ordering, is observed at low temperatures (2 - 19 K). 
At and above Tn long-range IC order predominates and a 
third order reflection is also observed as shown for T=20.99 
K. Above T ~ 21.7 K, broader peaks associated with the IC 
are observed up to « 36 K. B) Longitudinal scans, close to 
the origin, along the (0,K,0) direction show a peak compatible 
in position with the IC peaks observed at (0,1±(5,0). This 
peak, due to a ferromagnetic component of the modulated 
spin structure, is consistent with Model-I as described in the 
text (scans shown were obtained after subtraction of similar 
scans taken at 40K; scans are shifted in intensity for clarity). 

surate order persists in a narrow temperature range (~ 
0.9 K) above the C-IC transition. In this temperature 
range, third order reflections of the modulated IC struc- 
ture (see Fig. EI A), T = 20.99 K) and a peak close to the 
origin, characteristic of ferromagnetic modulations (Fig. 
|2lB)), are observed. No evidence for modulations along 
any other principal direction were observed. 

Figure O shows transverse scans at the IC peak 
(0,1±Q,0) with strong diffuse scattering below Tn, sig- 
naling the onset of the incommensurate phases at ele- 
vated temperatures. The onset for this diffuse scattering 
at T « 19K correlates with the maximum observed in 
the the ME-efFect, shown in Fig. [2 The incommensu- 
rate peaks are observed up to Tco ~ 36 K; however, in 
the temperature range 20.8 - to 21.70 K, they consist of 
two superimposed peaks, one resolution limited and the 
other diffuse. This leads us to conclude that LiNiP04 
undergoes two transitions, one from the short-range IC 
order to the long-range IC structure at T/c = 21.69, and 
a second at Tn from IC-LRO to AF as observed in the 
ME effect measurements. 

To account for the observations, two magnetic mod- 
els were considered in which each spin is rotated either 
about the b-axis or the a-axis by an angle a with re- 
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spect to its nearest neighbor (Model I or Model II, re- 
spectively). The angle of the j'th spin with respect to 
the spin at some arbitrary origin is given by aj — Q ■ Vj , 
where Q is a vector along the (010) direction as the 
AF propagation vector. The magnetic moment in the 
plane varies as follows, S = ^(sin Q ■ r,Q, e*'^ '' cos Q ■ r) 
for Model I and S = /i(0, sinQ • r, e''^ '' cos Q • r) for 
Model II. Whereas Model I predicts ferromagnetic mod- 
ulations with peaks near the origin at q = (0, ±(3,0), 
Model II does not. To determine the suitable model, 
scans close to the origin (at small q's) along all principal 
directions were conducted of which only the (0,K,0) scan 
gave evidence to the IC structure. Figure I2B) shows 
background-subtracted scans along the (0,K,0) direction 
(background was measured at T = 40 K) with a peak ob- 
served at the intermediate long-range IC phase (between 
T = 20.8 K and T = 20.7 K ). The intensities of the satel- 
lite peaks close to the AF propagation vector (see Fig. 
121) are about a hundred times stronger than the peaks 
near the origin at (0, ±(5,0), consistent with Model I as 
shown below. The intensity ratio I{t^p ± Q)/I{±Q) of 
these peaks can be estimated from the structure factors 
of the two reflections, Fm ^ J2j e^'^'^'Q x {S x q) with 
q = T^p± nQ; n = 1,3, yielding 



/(t,,±Q) \FiT,,±Q) 



H±Q) 



\Fi±Q)\ 



(1) 



where b is the lattice spacing along the direction of the 
modulation. Equation (^) shows the intensity of the peak 
near the origin vanishes as Q gets smaller, i.e., as the 
temperature is lowered towards the C-IC transition, giv- 
ing rise to a maximum in peak intensity as qualitatively 
shown in Fig. |21 

The IC magnetic structure occurs as an intermedi- 
ate phase between two high symmetry phases. As the 
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FIG. 3: A) Transverse scans along L B) and along H at 
roughly the position of the IC peak. In the temperature range 
20.8 - to 21.7 K (i.e., from Tn = Tc-ic to Tic) a broad dif- 
fuse peak is superimposed on a resolution limited peak (both 
Lorentzians) . Short range IC phase persists up to Too ~ 36 
K. 
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FIG. 4: A)Temperature-dependence of the order parameter 
in LiNiP04 as measured on the (0,1,0) magnetic Bragg re- 
flection. Intermediate phases between the paramagnetic (at 
temperatures higher than ~ 36 K) and the AF phases are in- 
dicated. In the IC region two phases are identified, one with 
long-range and one with short-range order. Diffuse scatter- 
ing at the nominal position persists up to Tco « 36 K (a 
crossover temperature). The transition from short- to long- 
range IC order occurs at Tic = 21. 7K, and the C-IC tran- 
sition at Tc-ic = Tn = 20.8 K. Illustration of the ground 
state magnetic ordering with a projection of the relevant ions 
on the b-c plane is also shown. B) Temperature dependence 
of the IC wave-vector. C) Simplified ground state (shaded 
area) beside the model of the IC structure within one plane. 
Here, each spin is rotated at a finite angle with respect to 
a neighboring spin about the a-axis (Model II in the text). 
Model I is similar to the one shown above, except the spins 
are rotated about the b-axis, giving rise to a FM modulated 
structure with a detectable peak at the origin (small angles). 



temperature is lowered from the paramagnetic phase, an 
onset for the IC occurs at Tco ~ 36 K, with a grad- 
ual increase in the wavelength of the modulation and 
the correlation length. At Tic = 21.7 K, the coher- 
ence length diverges and higher and new harmonics ap- 
pear until the modulation wavelength coincides with the 
high symmetry AF phase at Tn = 20.8 K (See Fig. 
2J). The origin of the IC structure may be induced by 
subtle charge distortions, or it could be innate to the 
spin Hamiltonian. Since no evidence for a structural in- 
commensurability was found, we hypothesize that the 
IC phases (long- and short-range) originate from spin- 
dimensionality crossover, i.e., from a continuous Heisen- 
berg (or XY) type model to an Ising-model. The ground 
state of the system is AF with no modulations, indi- 
cating a spin Hamiltonian that does not include strong 
terms that invoke an IC ground state (for instance, a 
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Dzyaloshinsky-Moriya term), as recently suggested by 
Kharchenko et al.^^. Thus, the occurrence of intermedi- 
ate magnetic incommensurate phases in LiNiP04 has all 
the characteristics of typical structural IC phases, which 
also manifest soft-phonon modes jSj- In preliminary spin 
wave studies of LiNiP04, a temperature-dependent gap 
(1.5-2 meV), and a minimum in the spin- wave dispersion 
curve was observed 26J . Theoretical predictions suggest 
that frozen magnons are possible in a 2D system with 
random-distribution couplings p^. as recently suggested 
for these systems 0,^3- 

In summary, our observations of a C-IC first-order 
phase transition, and the short-range order below Tn 
explain the abrupt jump and the maximum in the ME- 
effect of LiNiP04 (as observed in Fig. Our results 
however do not explain the anomalous temperature de- 
pendence of the ME coefhcient which may be due to ME 
contributions of the spontaneous toroidal moment, which 
is magnetoelectric in naturej28j. Recently, the anomalous 



temperature-dependence of the ME coefficient of weakly 
ferromagnetic/ferroelectric/ferroelastic boracites has, for 
the first time, been explained by considering the con- 
tribution of the symmetry-allowed spontaneous toroidal 
moment 29, 3,CL i31i] . Neutron scattering studies under 
applied magnetic field may shed light on the anomalous 
behavior of the ME-effect. We hope that our findings 
will stimulate theoretical microscopic studies on the na- 
ture of this IC phase observed in LiNiP04, but not in the 
isostructural LiCoP04. 
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